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Chapter 1 


Astrophysical and Cosmological Searches for 
Lorentz Invariance Violation 


Shantanu Desai 


Abstract Lorentz invariance is one of the fundamental tenets of Special Relativ- 
ity, and has been extensively tested with laboratory and astrophysical observations. 
However, many quantum gravity models and theories beyond the Standard Model 
of Particle Physics predict a violation of Lorentz invariance at energies close to 
Planck scale. This article reviews observational and experimental tests of Lorentz 
invariance violation (LIV) with photons, neutrinos and gravitational waves. Most 
astrophysical tests of LIV using photons are based on searching for a correlation 
of the spectral lag data with redshift and energy. These have been primarily car- 
ried out using compact objects such as pulsars, Active Galactic Nuclei (AGN), and 
Gamma-ray bursts (GRB). There have also been some claims for LIV from some 
of these spectral lag observations with GRBs, which however are in conflict with 
the most stringent limits obtained from other LIV searches. Searches have also been 
carried out using polarization measurements from GRBs and AGNs. For neutrinos, 
tests have been made using both astrophysical observations at MeV energies (from 
SN 1987A) as well as in the TeV-PeV energy range based on IceCube observations, 
atmospheric neutrinos, and long-baseline neutrino oscillation experiments. Cosmo- 
logical tests of LIV entail looking for a constancy of the speed of light as a function 
of redshift using multiple observational probes, as well as looking for birefringence 
in Cosmic Microwave background observations. This article will review all of these 
aforementioned observational tests of LIV, including results which are in conflict 
with each other. 


1.1 Introduction 


In special relativity, the speed of light, c, is a Lorentz invariant quantity. However, 
this ansatz is not valid in many theories beyond the Standard Model of Particle 
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Physics and also various quantum gravity and string theory models In these mod- 
els, Lorentz invariance is expected to be broken at ultra-high energies close to the 
Planck scale (Ep) ~ 10!° GeV), and the speed of light is consequently a function of 
photon energy [1]. Alternately, one could think of the vacuum refractive index as 
been different from unity in such models [2]. The energy-dependent speed of light 
can then be written as [3]: 
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where s+ = +1 denotes the sign of the Lorentz Invariance violation (LIV), corre- 
sponding to subluminal (sj = +1) or superluminal (s+ = —1); Egg denotes the 


quantum gravity scale where LIV effects kick in, and n is a model-dependent term 
and is usually equal to one or two, corresponding to linear or quadratic LIV, respec- 
tively. The correspondence between n and different theories of LIV is discussed in 
[4]. 

A plethora of searches for LIV have been carried out using photons, neutrinos, 
and gravitational waves. The astrophysical sources used for searches of LIV with 
photons include pulsars, Active Galactic Nuclei (AGN) and Gamma-Ray Bursts 
(GRBs). Most of these searches using AGNs and GRBs using spectral lags, which 
can be defined as the difference in the arrival times between high energy and low 
energy photons, and is positive if the high energy photons arrive earlier than the 
low energy ones [5]. For a source at cosmological distances, Equation 1.1 leads to a 
first-order differential time delay for signals as a function of energy [6]! : 
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where h(z) = = is the dimensionless Hubble parameter as a function of redshift. 


For the current standard ACDM model [7], A(z) = V 2m (1 +z’)? + Q4. The ex- 
pression for A(z) in other cosmological models can be found in [8]. 

The observed spectral time lag (At;5,) corresponding to an energy interval AE 
can be modeled as a sum of two independent delays: 


b. 
Atobs = fj — tn = Athy’ + Atry (1.3) 
where Afobs is the observed spectral lag; tn and t; represent the arrival times of 
photons of energies Ej, E; where Ej > Eu AUS is the intrinsic time delay due to 
astrophysical emission in the observer frame, and Atzzy due to LIV. 


The intrinsic time delay in the observer frame is given by 


Ato = (1 +2)Atin, (1.4) 


! All literature prior to [6] had (1 +z) missing in the numerator of Eq. 1.2 
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where the (1 +z) term accounts for the cosmological time dilation. The time delay 
due to the linear and quadratic LIV models can be obtained from = 1 and n = 2, 
respectively of the following equation [6], 


l+n E} -Er l (+z), 
Atiy = —S+ L —— — dz 1.5 
Ee ( 2H, | Ebon J Jo hE) v 


where Egg,» is the Lorentz-violating or quantum gravity scale, above which Lorentz 
violation becomes the dominant contribution, and Ep, E; represent the upper and 
lower energy intervals for photons arriving at times t, and t;, respectively. In Eq. 1.5, 
n= l and n = 2 denote linear and quadratic LIV models, respectively. We should 
point that the time delay in Eq. 1.5 corresponds to one explicit model of LIV [6], 
which assumes that the spacetime translations are not affected in any LIV model. AI- 
ternative LIV and doubly-special relativity models have also been proposed where 
quantum gravity affects spacetime translations [9] Therefore, for such models, 
Eq. 1.5 would not be valid. Here, we shall mostly discuss results which use the for- 
mulation in [6]. We also note that all results prior to the work in [6] had a missing 
(1 +z) term in the integrand of Eq. 1.5. Results on constraints on Lorentz violating 
standard model extension (SME) [10] from astrophysical observations (circa 2008) 
can also be found in [11]. We also note that many works report limit on photon 
mass obtained from dispersion relations while searching for LIV. Other techniques 
include looking for signatures of photon splitting and photon decay. 

This chapter will collate all the observational results on LIV using the afore- 
mentioned spectral lag method as well as various other astrophysical objects and 
observations such as Pulsars polarization observations (from AGNs, GRBs, CMB), 
Cosmology, neutrinos, and also a few terrestrial experiments. We should also men- 
tion that many other reviews focusing on observational tests of LIV have also been 
written over the years. A non-exhaustive list includes [12, 13, 14, 15, 16, 17, 18, 19, 
20, 21, 22, 23, 24, 25] (see also references therein). However, since this is a rapidly 
evolving field, this review will attempt to cover a few topics not mentioned in some 
of the aforementioned reviews. Corresponding reviews emphasizing the theoretical 
aspects of LIV can be found in [26, 27, 28, 29, 30, 31, 32, 20, 23, 9, 4] and will 
not be discussed here. A review of experimental tests of special relativity, which 
is the theoretical foundation behind Lorentz invariance can be found in [33, 34]. A 
compilation of limits on photon mass compiled using pulsars, fast radio bursts, solar 
wind, solar system planets, etc can be found in Particle Data Group [35] as well as 
[20] (See Table 3 therein) and are not reviewed here. Limits on LIV from gravita- 
tional Cherenkov radiation are discussed m [36]. This review will also not discuss 
the closely related topic of tests of Weak equivalence principle using astrophysical 
observations (for eg. [37]) or tests of variation of fundamental constants [38, 39]. 

The outline of this chapter is as follows. We shall review LIV searches using 
pulsars in Sect. 1.2; Active Galactic Nuclei in Sect. 1.3; neutrinos in Sect. 1.4; Cos- 
mological observations in Sect. 1.5; gravitational waves in Sect. 1.6. LIV searches 
with GRBs (which occupy bulk of this chapter) are discussed in Sect. 1.7. Searches 
with polarization measurements are reviewed in Sect. 1.8. We conclude in Sect. 1.9. 
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1.2 Pulsars and other Galactic sources 


Pulsars are rotating neutron stars, which emit pulsed radio emissions with periods 
from milliseconds to a few seconds and magnetic fields ranging from 10? to 10!4 
G [40, 41]. Pulsars have proven to be wonderful laboratories for tests and applica- 
tions of a whole range of topics in Physics and Astronomy [42, 43]. 

The first test of LIV using pulsars was carried out using EGRET observations of 
the Crab pulsar from 70 MeV to 2 GeV [44]. This work looked an energy dependent 
pulsar arrival time. No statistically significant variation of the pulse arrival time with 
energy was found. Subsequently, a 95% lower limit on the energy scale of quantum 
gravity Egg of 1.8 x 10!5 GeV was set [44]. A similar search was subsequently car- 
ried out from the Crab pulsar in the energy range 600 GeV—1.2 TeV by the MAGIC 
collaboration [45]. Since no such energy-dependent variation was found, upper lim- 
its were set for both superluminal (s+ = —1) and subluminal cases (s+ = 1). These 
limits are Egg > (4.5 — 5.5) x 10!7 GeV for linear LIV and Egg > (5.3— 5.9) x 1019 
GeV for quadratic LIV. Recently, sub-PeV observations of the Crab nebula from the 
Tibet [46] and HAWC [47] collaborations have also been used to set lower limits 
on LIV scale by looking for photon decay and photon splitting [48]. Based on the 
absence of photon decay and splitting, lower limits of Ezzy > 4.1 x 10!4^ GeV and 
Ep > 1.9 x 10? GeV, respectively were set [48]. The HAWC collaboration also 
did a similar analysis of the Crab pulsar and three other pulsars to set limits on LIV 
with lower limits of (100—200) E,,; [49]. The Crab pulsar has also been observed 
at over 1 PeV by the LHAASO-KN2a collaboration with the most energetic pho- 
ton having an energy of 1.12 +0.09 PeV [50]. Due to LIV, a free electron can emit 
Cherenkov radiation in vacuum. Assuming the inverse Compton spectrum of the 
Crab nebula exceeds up to PeV energies, a constraint on LIV was obtained based 
on the stability of > 1 PeV electrons against electron Cherenkov radiation. Thus, a 
lower limit on Egg > 1025 GeV was obtained [51]. 

In 2021, the LHAASO collaboration detected 12 gamma-ray sources with en- 
ergies from 100 TeV up to 1.4 PeV with a statistical significance of > 70 [52]. 
Except for the Crab pulsar, these sources have not been confirmed from previously 
known point sources. A search for superluminal LIV was done using this dataset 
by looking for single photon pair production as well as single photon splitting into 
multiple photons [53]. Both these processes lead to a cut-off in the spectrum. The 
LIV parameters were obtained using a forward modelling procedure. No spectral 
cutoff due to LIV was seen and therefore a lower limit on LIV was set. The 9596 
c.l. lower limit on first order LIV is Egg > 10 E, [53]. The limits on second order 
LIV is given by Egg > Dier pl [53]. An independent search for superluminal LIV 
of this data was also carried out [54], which obtained approximately similar lower 
limits as in [53]. Furthermore, it was also pointed out that the observation of the 
highest energy photon (1.4 PeV) could be compatible with subluminal LIV with a 
scale of 3.6 x 10!7 GeV [54]. 

A search for LIV in the context of quantum electrodynamics was also done using 
the Crab pulsar [55]. In this setup, the maximum vacuum speed for an electron (ce) 
is different from c and can be written as: 
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Ce = cy(1 +6) (1.6) 


where |6| << 1. For 6 < 0, photon decays to an electron-positron pair for pho- 
tons with energies Emax = me./2/|6|. From the observations of Crab pulsar at 50 
TeV [56], we get 6 > —2 x 10716 [55]. If 5 > 0, electrons become superluminal for 
energies E > 0.5Emax and will emit radiation at all frequencies due to Cherenkov ra- 
dation. Based on the observations of electrons in primary cosmic-ray with energies 
up-to 1 TeV [57], one gets Emax > 2 TeV implying 5 < 1.3 x 107" [55]. 


1.3 Active galactic nuclei 


All tests of LIV using AGNs have been carried out using observations of blazars 
at TeV gamma rays with atmospheric Cherenkov telescopes. Blazars are AGNs 
with jets pointed directly at us and with Lorentz factors of about 10 [58] The key 
idea is to look for correlation of spectral lag with energy to constrain Egg using 
Eq. 1.5. The first test of LIV using Active galactic nuclei (AGN) observations was 
made based on the TeV gamma-ray flare from Mrk 421 on 15 May 1996 detected 
by the Whipple atmospheric Cherenkov telescope [59]. By comparing the arrival 
times of photons with energies < 1 TeV and > 1 TeV, a 95% c.l. lower limit of 
Egg > 0.4 x 10!7 GeV was set. Subsequently, the MAGIC collaboration analyzed 
two TeV flares of Markarian 501 (located at z — 0.034) between May and July 
2005 in the energy range from 0.15 - 10 TeV. Based on the time delay, a lower 
limit of Egg > 0.21 x 1015 GeV and > 0.26 x 10!! GeV for linear and quadratic 
LIV, respectively was obtained [60]. Note that these limits are valid for both sub- 
luminal and superluminal propagation. The HESS Collaboration then did a similar 
analysis looking for photon time delays using the TeV gamma-ray flare observed 
from PKS 2155-304 on July 28, 2006. Two different methods based on wavelet and 
the method of cross-correlation were used to estimate the time delays between the 
light curves of different energies. This method then yielded 95% c.l. lower limits 
of Egg > (5.2 — 7.2) x 107 GeV for linear LIV and Egg > 1.4 x 10? GeV for 
quadratic LIV, which are valid for both superluminal and subluminal effects [61]. 
This result was subsequently superseded by a refined analysis technique consisting 
of an event by event likelihood fit as applied to the flare data and new lower limits of 
Egg > 2.1 x 1015 GeV and Egg > 6.4 x 10!° GeV at 95% c.l. were obtained [62]. 

The HESS collaboration also used observations of a VHE flare with E > 100 
GeV from PG 1553+113 (located at a redshift z = 0.49 + 0.05) observed on the 
night of April 26 and 27, 2012 to constrain LIV [63]. No energy dependent delay in 
the photon arrival time was seen. Consequently, a lower limit of Egg > 4.11 x io” 
GeV and Egg > 2.1 x 10!° GeV was set for linear and quadratic LIV, respectively, 
applicable for both superluminal and subluminal propagation [63]. 

An independent analysis of the above flare has also been carried out [64]. This 
work obtained lower limits of Egg > 4.9 x 1016 GeV, 1.2 x 10!” GeV, and 2.6 x 1018 
GeV for Mrk 421, Mrk 501, and PKS 2155-304, respectively for linear LIV [64]. 
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For a quadratic model, the corresponding limits obtained were Egg > 1.5 x 1010 
GeV and 9.1 x 10!° GeV for Markarian 421 and PKS 2155-304, respectively [64]. 

The HESS collaboration detected another TeV flare from Mrk 501 in June 2014 
with multi-TeV variability on time scales of minutes with its energy spectrum spec- 
trum extending up to 20 TeV [65]. This data was used to search for evidence of LIV 
for both s+ = +1 as well as s+ = — 1. The best-fit lower limit for linear LIV model is 
given by Egg > 2.6, 3.6(x 10^ GeV for s+ = 1 and s+ = — 1, respectively [65]. The 
corresponding limits for quadratic LIV model are given by Egg > 8.5,7.3(x 1019) 
GeV for s+ = 1 and s+ = — 1 respectively [65]. 

Some of the earlier limits on LIV obtained using Mrk 421, Mrk 501, and PKS 
2155-304 were disputed in [66]. This work showed that with a careful re-analysis 
of the light curves, the earlier upper limits should be revised and are consistent 
with light speed variation as a function of energy, which were hinted at from GRB 
observations (to be discussed in the subsequent section) [66]. More details on these 
methods and comparison with the previous limits on LIV for the same object can be 
found in [66]. 

In addition to the aforementioned searches for LIV from single sources, lim- 
its have also been set using a stacked catalog of blazars based on the distortion in 
the observed spectra, which gets accumulated during propagation. Two examples of 
such an approach include Ref. [67] and Ref. [68]. The key idea behind this is that, 
in the presence of LIV, the threshold for pair production gets lowered [6], which 
affects the optical depth of > 10 TeV photons from extragalactic sources. In the 
first such study, 86 spectra from 38 AGNS were used [67]. A joint fit to the ab- 
sorption signature seen in gamma-rays using models of extragalactic background 
light (EBL) along with Egg was done. The first work which applied this technique 
obtained a 95% c.l. lower limit on Egg given by Egg > 0.6E,,; for linear LIV and 
s+ = 1 [67]. This method of stacking spectra from multiple AGNs was then im- 
proved by selecting the optimum bins to account for the uncertainty in extragalactic 
background light [68]. They applied this technique to 18 measured spectra from 
six AGNs detected at TeV energies from HESS, HEGRA, VERITAS, and TAC- 
TIC for sub-luminal propagation (corresponding to s+ = 1). The 95% c.l. lower 
limit on Egg for linear and quadratic LIV are given by Egg > 12.8 x 107? GeV 
and Egg > 2.38 x 10!2 GeV, respectively [68]. These constitute some of the most 
stringent limits on LIV to date. 

A search for LIV within QED was also carried out with observations of Mrk 
501 (along the same lines as described in Sect. 1.2 for Crab pulsar), based on the 
fact that there is no decrease in the inferred optical depth based on the observed 
spectrum [55]. This work obtained a limit on ô given by ô < 2 x 10-16 [55]. 


1.4 Neutrinos 


Limits on LIV from neutrinos have been obtained from both neutrino oscillation 
experiments based on atmospheric neutrinos as well as of astrophysical neutrino 
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observations, and also using direct tests of speed of neutrinos using long baseline 
neutrino oscillation experiments. We briefly summarize some of these key results. 
Other reviews of tests of Lorentz invariance and CPT violation using neutrinos can 
be found in [69, 70]. 

The year 1998 was a watershed moment in the history of Physics when evidence 
for neutrino oscillations was found using the Super-Kamiokande experiment from 
atmospheric neutrinos [71], and subsequently re-affirmed with more statistics and 
using all the combined atmospheric neutrino samples [72, 73, 74]. After the first 
evidence for atmospheric neutrino oscillations [71], LIV was also considered as 
a possible mechanism for these neutrino oscillations [75]. For oscillations due to 
LIV, the oscillation probability scales as L x E as opposed to L/E for atmospheric 
neutrino oscillations [76]. However the oscillation probability using the 1998 data 
scales with LE as (LE) 9?*9^ [76] indicating that the oscillations seen in Super- 
K cannot be due to LIV. In 2015, the Super-K collaboration carried out a dedicated 
search for LIV using the full neutrino dataset from 100 MeV to over 100 TeV and did 
not find any evidence for the same [77]. Consequently, constraints on coefficients of 
Lorentz violating SME were set. A similar search was also done using the IceCube 
neutrino dataset between 100 GeV and 1 PeV [78]. 

In 2011, there was a flurry of excitement when the Opera experiment claimed 
evidence for superluminal neutrinos. However about a year later, new sources of 
errors were found in the original measurements and the refined analysis showed that 
there is no evidence for supeluminal neutrinos [79]. Other long baseline neutrino 
oscillation experiments such as ICARUS [80] and MINOS [81] have also found 
no evidence for superluminal neutrinos. More details on the searches and potential 
implications for superluminal particles can be found in a recent review [82]. 

The first (and to-date the only) observation of astrophysical MeV neutrinos from 
SN 1987A from the LMC at a distance of 50 kpc opened up the era of multi- 
messenger astronomy with neutrinos [83, 84]. It was shown from the simultaneous 
arrival of photons and neutrinos after a journey of 160,000 years, that the velocity 
of neutrinos is comparable to the speed of light with an accuracy of 8 (10^?) [85]. 
Similar constraints on the neutrino speed were obtained using observations of PeV 
neutrinos from TXS 0506+056 by IceCube with difference between neutrino speed 
and speed of light equal to within CU 19-15) [86] (See also [87]), which is six orders 
of magnitude more stringent than SN 1987A. From the same event a constraint of 
Eog > 3 x 10!6 GeV was obtained [88, 89]. 

Searches for LIV using astrophysical neutrinos have also been done using neutri- 
nos with energies in the TeV-PeV range. Following the first detection of two 1 PeV 
cascade events by the IceCube collaboration [90], a lower limit of Egg > 10°E pl 
and Egg > ed? pi Was obtained for linear and quadratic LIV, respectively assum- 
ing both the events are extragalactic [91]. There have also been a class of searches 
which have set limits on LIV by looking for association between GRBs and IceCube 
neutrinos [92, 93, 94, 95, 96]. We provide some highlights from these works. More 
details can be found in the original references. Ref. [94] found 13 IceCube neutri- 
nos with energies between 60 TeV and 2 PeV, which are associated with GRBs. 
They showed that the arrival time of neutrino events is consistent with an energy- 
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dependent speed of light variation indicating positive evidence for LIV with an en- 
ergy scale of Eu = (6.4 x 1.5) x 10!” GeV. Ref. [96] looked for a correlation 
between the high energy starting events in 7.5 years of IceCube data with more than 
7000 GRBs. This work found six possible GRB-neutrino pairs. They used a variant 
of Eq. 1.5 as follows [96]: 


Apa, | Sve SgBy is (14- z)dz (1.7) 
Ho Eog 0 Gell Lt QA 


where z is the GRB redshift, Ey, E, are the neutrino and photon energies respec- 
tively; sy and s, are the sign factors for the photon and neutrino, respectively. 

Based on this, a lower limit of Egg > 8 x 10! GeV was set. One issue with 
the above (and similar other results [92, 93]) is that no thorough statistical analysis 
to assess the significance of the association between GRBs and neutrinos has been 
done. It is entirely likely that the neutrinos could be due to atmospheric neutrino or 
diffuse astrophysical neutrino background (not associated with GRBs). The analysis 
done by the IceCube collaboration does not find any evidence for neutrinos spatially 
or temporally correlated with GRBs [97, 98]. Nevertheless, as IceCube is continu- 
ing to take data and with future upgrades such as IceCube-Gen2 [99] as well as 
PINGU [100], it may be possible to detect ultra-high energy neutrinos from GRBs, 
which would shed additional light on LIV. 


1.5 Cosmological observations 


The standard (ACDM) Cosmological model consisting of 7096 dark energy and 
25% dark matter agrees well with large scale structure observations [7]. Recently a 
number of observational tensions have been found in this standard model such as Ho 
tension, Og tension etc [101]. The edifice of the standard model is General Relativity 
which subsumes Special Relativity. Therefore, Lorentz Invariance and the constancy 
of speed of light is inbuilt in the Standard cosmological model. Although no one has 
thoroughly explored the implications of LIV based theories to resolve or alleviate 
some of the cosmological tensions, many tests of alternatives of ACDM, which posit 
a varying speed of light as a function of redshift (c = c(z)) have been carried out. 
Some of these theories can solve the cosmic coincidence problem. One such class of 
models consist of Varying speed of light theories [102, 103]. Many tests of varying 
speed of light theories have been carried out using cosmological observations using 
BAO, galaxy clusters, strongly lensed sources along with Type la SN data. The key 
idea behind testing the variability of speed of light is to test for variation of D,(z): 


1 z cdz 
Da(z) = rah Hz) (1.8) 


where c is the speed of light. Eq. 1.8 can be generalized by replacing c with a vary- 
ing c(z). A cosmological measurement of the variability of speed of light has been 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 9 


carried out using the Pantheon Type 1a Supernovae sample, Baryon Acoustic oscil- 
lation measurements and cosmic chronometer observations [104]. This work found 
a constant speed of light at two different redshifts c = (3.2 +0.16) x 10° km/sec 
at z 1.58 and c = (2.67 +0.14) x 10° km/sec at z ~ 1.58 at z © 1.36. A similar 
study with galaxy clusters done by combining the measurements of galaxy cluster 
gas mass fractions, H(z) from cosmic chronometers, and Type-Ia Supernova also 
found no variation of speed of light with redshift [105]. Another test for the vari- 
ation of the speed of light was done by combining data from four strongly lensed 
systems (SLACS, BELLS, SL2S and LSD) along with Type 1a SNe [106]. No de- 
viation in the speed of light from c as found as a function of redshift. Various other 
astrophysical and cosmological tests of varying speed of light theories are reviewed 
in [107]. 

Another consequence of LIV is the rotation of linear polarization also known 
as birefringence, where the photon dispersion relation depends on the polarization. 
This causes a rotation of the Cosmic Microwave background (CMB) polarization 
and induces parity odd cross-correlations [108] The first such search was done using 
the rotation measure limits from 5 years of WMAP data [109], and a limit on the 
effective photon mass was found with my < 3.8 x 10-?ev [110]. Recently, a 2.40 
evidence for non-zero birefringence was found using Planck public release 3, with 
the birefringence angle (f) given by B = 0.35 +0.14 [111]. A refined value of 
P — 0.30 +0.11 was obtained using Planck public release 4, with the significance 
getting enhanced to 3o [112]. Further implications of this result and prospects for 
confirmation with upcoming experiments are reviewed in [113]. 

Other tests with CMB data include signatures of specific LIV theories such as 
effect of Einstein-Aether theory on the CMB anisotropy spectrum [114] and bire- 
fringence due to Lorentz violating Standard model extension [115]. 


1.6 Gravitational Waves 


The year 2016 ushered in the era of gravitational wave (GW) astronomy with the 
detection of binary black hole merger GW150914 by the LIGO-VIRGO collabora- 
tion [116]. To-date the LIGO-VIRGO collaboration has detected GWs from about 
100 binary black hole, binary neutron star or black hole/neutron star binaries. Tests 
of LIV with gravitational waves entail a direct test of speed of gravitational waves. 
A direct test of speed of gravitational waves has been done from the first binary 
black hole merger which showed that c,,, < 1.7c [117]. The simultaneous detection 
of GWs and photons from GW 150914 also enabled a direct test of relative vac- 


uum velocity between GWs and photons —3 x 10715 < ( = ) < 7 x 10716 [118]. 


C 


Other indirect tests based on Shapiro delay of GWs can be found in [37, 119]. 
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1.7 Gamma-Ray Bursts 


The vast majority of tests of LIV have been done using GRBs for more than two 
decades starting with [3], which did the first feasibility study demonstrating the 
power of using GRBs for probing LIV violation. Doing a thorough comprehensive 
summary of over two decades of literature is beyond the scope of this short review. 
So we shall mainly report highlights each of the works which used GRBs from 
LIV. More details can be found in the original references or in some of the other 
observational reviews on LIV mentioned in the Introduction. 

GRBs are single-shot explosive events, which have been observed over ten 
decades in energies from keV to 100 TeV range [120]. They are located at cos- 
mological distances, although a distinct time-dilation signature in the light curves 
has not yet been demonstrated [121]. GRBs are traditionally divided into two cate- 
gories based on their durations, with long (short) GRBs lasting more (less) than two 
seconds [122]. Long GRBs are usually associated with core-collapse SN [123] and 
short GRBs with neutron star mergers [124]. There are however many exceptions 
to the aforementioned dichotomy, and arguments for additional GRB sub-classes 
have also been made [125, 126] (and references therein). Similar to AGNs, most 
searches of LIV with GRBs have beend one using spectral lags. Besides spectral 
lags, a few tests have also been done using polarized data from GRBs. In this sec- 
tion we discuss results from spectral lag searches starting with fixed observer frame 
energies and then using fixed energy intervals in the source frame. In Sect. 1.8, we 
shall discuss LIV searches with GRBs using polarization measurements. We note 
that a tabular summary of results for LIV using GRB spectral lags (upto 2021) has 
also been collated in Table 1 of [21]. 


1.7.1 Searches in fixed observer frame energies 


The first study for LIV using GRB spectral lags was carried out using GRB 930229 
based on the simultaneous arrival of a flare with a rise time of (220 + 30) u-sec 
that occurred simultaneously from 30 — 200 keV [127]. From this GRB, a limit on 
fractional variation of speed of light given by Ac/c < 6.3 x 10?! was set [127]. 
Subsequently, data from 6 GRBs with measured redshifts from the BATSE as well 
as OSSE detector was analyzed to look for energy dependent speed of light by look- 
ing for a correlation in the spectral lag with redshift and also stochastic variation 
in photon velocities at the same energies [128]. No such variation was found and 
subsequently a limit of Egg > 107? GeV was set [128]. A refined analysis using 
wavelets increased the lower limit to Egg > 6.9 x 10!5 GeV (for subluminal propa- 
gation) [129]. The first test of LIV using GRB spectral lags in the MeV energy range 
was carried out using three different energy bins from 0.2-17 MeV using data from 
RHESSI satellite [130]. No time lags was found. The limit on Egg was obtained 
from Eq 1.2, assuming the source is at redshift of z ~ 0.3, although this was not a 
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direct redshift estimate. This work obtained a limit of Egg > 1.8 x 10!” GeV and 
5.5 x 107 GeV for linear and quadratic LIV, respectively [130]. 

The first such systematic study by stacking a GRB sample was carried out in 
[2, 131], who considered a sample of 35 GRBs in the redshift range z = 0.168 — 4.3 
from HETE, BATSE, and Neil Gehrels SWIFT. This was the first work where the 
spectral lags were modeled as the sum of a constant intrinsic lag along with another 
contribution due to LIV (Eq. 1.3). This was also the first work where an extra (1 + 
z) term was used inside the integrand in Eq. 1.5 following the arguments in [6]. 
They also found a 40 evidence for the higher energy photons to arrive earlier than 
the lower energy ones, which at face value points to evidence for LIV with s+ = 
—] [2]. However, when an additional systematic offset was added to enforce the 
x°/DOF to be equal to one, the statistical significance for LIV reduced to about lo. 
Subsequently, a lower limit of Egg > (0.9 — 2.1) x 10/6 GeV was set at 95% c.1 [2]. 
Note that even though at face value, the results showed tentative hints for e = —1, 
the limits in [2] were obtained for s+ = 1. The same data was also used to search for 
LIV by combining it with BAO and Union2 Type 1a SN data, and within the aegis of 
three cosmological models, namelyA CDM, wCDM, and Chevallier-Polarski-Linder 
model [132]. The results were consistent with a zero slope indicating that there is 
no evidence for LIV [132]. A joint cosmological and LIV analysis of the same data 
along with JLA Type la SN data has also been carried out, which showed that the 
slope is also consistent with no LIV to within 20 [133]. 

Subsequently, a similar test for spectral lag variation with energu was done us- 
ing GRB 0512214, located at z = 0.5465, using light curves from Swift-BAT, and 
Konus-Wind [134]. For the Swift data, the spectral lags were estimated using four 
different energy bins 15-35 keV, 50-150 keV, 100-150 keV, and 300-350 keV. Based 
on the Swift time resolution of 4 msec, no time lags were detected between the 
peaks in the different energy bands and a limit on LIV energy scale was set with 
values given by (0.01 — 0.033)E,; for linear LIV and (1.2 — 3.6) x 107 PE, for 
quadratic LIV for s+ = 1 [134]. For the Konus-Wind data (having a time resolution 
of 2 msec), spectral lags were calculated based on the 18-70 keV, 70-300 keV, and 
300-1160 keV energy bands. No time delays between the different energy bands 
were detected, and subsequently limits on LIV were set whose values are given by 
(0.0012 — 0.0066)E,; for linear LIV and (1.1 — 5.1) x i0 E, for s+ = 1 [134]. 
The slight variation in the aforementioned limits changes depending on which pair 
of energy bands were used to calculate the spectral lags. In this work, more stringent 
limits were obtained from Konus-WIND data due to better time resolution. 

A subsequent search was then done using Integral data [135]. This work used un- 
binned maximum likelihood analysis of 11 GRBs with known redshifts to measure 
17 time lags. A fit to the time delay consisting of a constant intrinsic lag along with 
LIV induced term was done. However the best-fit Y?/dof was much greater than 1. 
After rescaling the error a lower limit of Egg > 3.2 x 10!! GeV was found at 9596 
c.l [135]. This limit is valid for both subluminal and superluminal propagation. 

Then a similar search was carried out using HETE-II data [136]. They considered 
15 GRBs in the redshift range z — 0.16 — 3.37 in the energy range of 6-400 keV. 
The spectral lags were calculated using de-noised light curves between 8-30 keV 
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and 60-350 keV. A search for time lag as a function of both redshift and energy 
difference was carried out. No significant quantum gravity signature was seen. A 
constant intrinsic time lag was considered and data fit to Eq. 1.3. A lower limit of 
Eog > 2 x 10/5 GeV was obtained at 95% c.l. for s+ = 1 [136]. 

About a month after the launch of the Fermi satellite, a bright GRB (GRB 
080916C) located at a redshift (located at z — 4.35 +0.15)was detected [137]. Based 
on the spectral time lag of 16 seconds between a photon of energy 13 GeV and low 
energy photons of keV energy, a lower limit of Egg > 1.2 x 10'° GeV was obtained 
for s+ = 1 [137]. Subsequently, a LIV search was carried out using GRB 090510 lo- 
cated at z = 0.903 + 0.003 by calculating the spectral lag between photons above 30 
GeV and below 1 MeV [138]. This gives a constraint of Egg > 7.72 x 10!? GeV and 
Egg > 7.26 x 10!° GeV for linear and quadratic LIV for s+ = 1, respectively [138]. 
A similar analysis was also done by the Fermi LAT/GBM collaboration and multi- 
ple lower limits of (10 — 100)Ep; at 9596 c.l. were obtained for linear LIV for both 
values of s+ (cf. Table 2 of [139]), where the limits change slightly depending on 
which energy bands were used to calculate the spectral lags [139]. This was the first 
work which obtained a limit on Egg greater than the Planck scale for superluminal 
and subluminal LIV. 

Another search for LIV was then done by considering 4 GRBs from Fermi-LAT 
spanning the redshift range from 0.9—4.13 [64]. Similar to [2]. they assumed that 
the spectral lag consists of a constant time lag along with a LIV-induced lag. This 
work recast Eq. 1.3 and Eq. 1.4 into a regression relation of the form (assuming 
$4 — —1) 


em = Atin + K/Epg (1.9) 
where K = Ain: X Epc- A separate fit was done for the three long GRBs as well 
as one short GRB in the sample. They showed that the spectral lag of the three 
long GRBs can be well fitted by the sum of a constant lag along with a linear LIV 
lag [64]. This was the first work using GRB spectral lags, which found evidence 
for LIV with 1o bounded intervals for Egg given by: Egg = (2.2 + 0.9) x 10" 
GeV and Egg = (5.4+0.2) x 10° for linear and quadratic LIV, respectively [64]. 
However they pointed out that if assumes that the observed lag has contributions 
from LIV, one obtains negative values for the intrinsic astrophysical lag, which is 
contrary to conventional wisdom for the astrophysical models, if low energy photons 
are produced by electrons and high energy ones by protons. Therefore, they urged 
that these hints for LIV should be treated with caution [64]. 

A somewhat novel search for LIV using Fermi-LAT GRBs was done by mea- 
suring the spacing between photon bunches in four GRBs, which were found to 
be much shorter (for GRB 090510A) than any other temporal feature in the en- 
ergy range from 1-30 GeV. [140]. Based on this the most stringent bound on 
Ac/c « 6.94 x 107?! was obtained for GRB 090510A [140]. 

Then, Vasileiou et al [141] considered 4 Fermi-LAT GRBs in the GeV energy 
range with measured redshifts. This work used three complementary techniques to 
quantify the amount of spectral dispersion in the data. Then using conservatives 
choices for the source intrinsic spectral evolution, constraints were obtained on LIV. 
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For s4 = 1, the most stringent limit was obtained for GRB 090510 with Egg > 
7.6E for linear and Egg > 1.3 x 10!! GeV for quadratic LIV [141]. For s4 = —1 
the corresponding best limit was Egg > 18E,; using GRB 080916C for linear and 
Egg > 6.7 x 10!° GeV using GRB 090510 for quadratic LIV [141]. The same group 
of authors then analyzed Fermi-LAT data from GRB 090510 to look for a stochastic 
spread in the speed of light around c [142]. No such variations were seen and 95% 
c.l. lower limits on Egg are give by Egg > 2.8E, [142]. 

The next work which confirmed the previous tentative evidence for LIV found in 
[64], was by [143], which followed the same procedure as in [64, 2]. They analyzed 
the spectral lag data for 8 GRBs detected by both Fermi-LAT and Fermi-GBM. 
They assumed a constant intrinsic lag in the source frame and fitted the spectral lag 
to Eq. 1.9. They showed that data from 5 long GRBs could be fit using a straight 
line with an intrinsic lag of ting = —12.1 + 1.7 secs and Egg = (3.05 +0.19) x io" 
GeV for linear LIV hypothesis. After taking into account the uncertainty in the 
cosmological parameters and energy resolution of LAT, this uncertainty in Egg is 
equal to 0.7 x 10!7 GeV [143] If we consider a quadratic LIV, the best-fit value of 
Egg = (6.8 + 1.0) x 10? GeV [143] with s+ = —1. A variant of this analysis using 
a subset of the GRBs analyzed in [143] along with a few additional ones confirmed 
this tentative evidence for LIV with Egg = (3.6 +0.26) x 107 GeV [144]. More 
confirming evidence was found using GRB 1605094 (located at z = 1.17) using 
photons in the energy range 1—50 GeV with Egg — 3.6 x 10!7 GeV [145]. Two 
new methods for characterizing the arrival times of low energy photons needed for 
calculation of spectral lags in Fermi-GBM and Fermi-LAT data were proposed in 
[146] and applied to data previously analyzed [144, 145]. Although there was a 
slight difference in the regression relation, the results in [146] were consistent with 
those in [144, 145] 

The first search using a sample consisting of only short GRBs detected by 
SWIFT/BAT and with redshifts in the range 0.36—2.2 was carried out in [147]. The 
rationale for selecting only short GRBs is that these GRBs have negligible spectral 
lag and smaller intrinsic dispersions compared to long GRBs. Therefore, choosing 
an exclusive sample of short GRBs would minimize any systematics related to the 
unknown nature of the intrinsic lag. The spectral lags were estimated between the 
50-100 keV and 150-200 keV energy intervals. This dataset of 15 GRBs was fit to 
a constant intrinsic lag and a LIV-induced lag (cf. Eq. 1.9). While fitting the data an 
extra scatter was added to account for the dispersion in the spectral lag. However, 
unlike previous works [64, 143, 145, 144], the slope of the regression relation in 
Eq. 1.9 was consistent with zero, implying no LIV. Assuming a linear model for 
LIV, a 95% c.l. lower limit on Egg was set given by Egg > 1.48 x 1016 GeV, valid 
for both values of s+ [147]. 

The first convincing case for a spectral lag transition from positive to negative 
lags, using multiple lags over three decades in energy range from the same GRB 
was demonstrated for GRB 160625B [148] located at z — 1.41. This burst had three 
different sub-bursts The spectral lags were calculated for 30 time lags starting from 
12-16 keV to 17-20 MeV. This was also the first work, which did not assume a 
constant intrinsic lag model, and instead used a phenomenological model for the 
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intrinsic time lags given by: 


TERORO I 


where a depicts the energy exponent and T depicts the time scale for the intrinsic 
time lag, Eo and E correspond to the lower and upper energy intervals. We note that 
if the energy interval corresponds to finite energy band, the midpoint of that band 
is usually used to ascribe an energy to it. This has been the norm in all literature 
on LIV searches using spectral lags. However, technically one should include the 
energy spectrum weighted values. This intrinsic model was empirically determined 
by modelling the single-pulse properties of about 50 GRBs [149]. This model has 
also been used for blazar flare modelling [150]. 

The analysis in [148] found that the spectral lag for this GRB shows a turnover 
at 8 MeV, with a transition from positive to negative lags. They argued that this 
transition could be a signature of LIV, which dominates at high energies. Their 
best-fit value for Egg was log(Egg/GeV) are 15.667055 and 7471000 for linear 
and quadratic LIV models, respectively assuming s+ = 1 [148]. In a subsequent 
work, the statistical significance for this spectral lag transition was then ascertained 
using both frequentist [151] and information theory techniques such as AIC and 
BIC [152, 153]. The frequentist test corresponds to a Z-score between 3.05 — 3.746, 
and AAIC/BIC > 10 for the quadratic LIV model, which point to decisive evi- 
dence [154] (see also [155]) However, despite finding a signature for LIV and ob- 
taining marginalized closed contour intervals for Egg, in the conclusions, only one- 
sided lower limits for Egg were reported in the conclusions, (Egg > 0.5 x 1016 GeV 
and Egg > 1.4x 10’ GeV, for linear and quadratic LIV, respectively for s+ = 1. Con- 


straints on the isotropic and anisotropic coefficients of the Lorentz violation SME 
using this GRB were also obtained [156]. However, subsequently it was pointed out 
that the spectral lags are correlated with the spectral evolution in GRB pulses [157]. 
This turnover in the spectral lag data could therefore be the consequence of a radiat- 
ing relativistic jet shell with a cutoff Band power law spectrum, similar to that seen 
in observations [158]. This model was also used to explain the spectral lag tran- 
sition in GRB 1905304, where however they used a three segment broken power 
law [159]. 

Subsequently, the spectral lag data for GRB 160625B was stacked together with 
the data for 35 GRBs from [2] to look for LIV corresponding to s+ = —1 [160]. Us- 
ing the same intrinsic model as [148] this stacking procedure provides a more robust 
estimate of the intrinsic time lag. The best-fit estimates obtained are log(Egg / GeV) 
are 14.523 ee and 8.79 + 0.0097 for linear and quadratic LIV models, respec- 
tively [160]. This was also the first work which used a model-independent estimate 
of the expansion history using Gaussian process regression [161, 162] based on H(z) 
measurements from cosmic chronometers [163]. 

Most recently, a similar spectral lag transition from positive to negative lags 
(similar to GRB 1606025B) was also detected in GRB 1901114C (detected by 
SWIFT) at about 0.7 MeV [164]. This work calculated the spectral lags in en- 


ergy intervals from 15-5000 keV and compared them to the lowest energy value 
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of 12.5 keV. However, the model for the intrinsic time lag asserted was opposite in 
sign compared to [148, 160]. The statistical significance of the spectral lag transi- 
tion was also estimated using Bayesian model comparison and found to be > 100, 
pointing to decisive evidence using Jeffrey’s scale. Similar to [148, 160], they ob- 
tained best-fit values for linear and quadratic models, with their 20 bounds given by 
log(Egg/GeV)= 14.49*012 and 6.00 + 0.06, respectively assuming s+ = 1. They 
showed using X?/DOF that all the LIV models provide a good fit to the LIV hypoth- 
esis [164]. 

The MAGIC collaboration however failed to find a similar evidence for an 
energy-dependent speed of light in the TeV gamma ray data for the same GRB. Us- 
ing conservative assumptions on spectral and temporal evolution, the MAGIC col- 
laboration set a 9596 c.l. lower limit on Egg of 0.55 x 10! GeV and 5.6 x 101? GeV 
for linear and quadratic models, respectively for superluminal propagation [165]. 
For subluminal propagation, the corresponding limits are 0.28 x 101? GeV and 
7.3 x 101? GeV for linear and quadratic models, respectively [165]. 

Therefore, the situation as of 2021 was that some of the above aforemen- 
tioned analyses obtained closed lo confidence intervals for Egg, which is less 
than the Planck energy scale (indicating evidence for LIV), although not with 
the same value of s+. Furtheremore, the estimated LIV energy scale estimated 
in [148, 160, 164] contradicts the previous the most stringent lower limits on 
LIV [139, 141, 142]. Still, given the tantalizing hints for LIV in some of the afore- 
mentioned works [2, 148, 160, 164], if the spectral lag data for all the GRBs can 
be described by a unified model, one would expect the net statistical significance 
of the LIV to get enhanced, if one stacks the spectral lag data from all the GRBs, 
and analyze them uniformly with the same model. Therefore, in order to resolve 
this imbroglio, the data from all the works which found these hints for LIV signa- 
tures in the spectral lag data [2, 148, 160, 164] were stacked together [166]. In all, 
there were a total of 91 spectral lag measurements from 37 GRBs, including 19 lags 
from GRB 190114C and 37 lags from GRB 1606025B. The aforementioned work 
did a Bayesian model comparison that the combined spectral lag data is a combina- 
tion of both an intrinsic and a LIV-induced lag, as opposed to only an intrinsic lag. 
The model assumed for the intrinsic astrophysical emission is the same as that pro- 
posed in [148, 160]. The priors used for œ and T were broad enough to encompass 
both positive and negative values. The expansion history was estimated in a non- 
parametric fashion using Gaussian process regression. When searching for LIV, no 
closed contour was obtained for linear LIV (see Fig. 2 of [166]). Only the quadratic 
LIV showed a closed contour. They found that the Bayes factor for both the LIV 
models was greater than 100 relative to the null hypothesis. The 68% c.l. lower limit 
for Egg was log[Egg(GeV )] > 16.08 for linear LIV, whereas for quadratic LIV the 
best-fit value obtained was log/Egg(GeV)| = 7.17 + 0.07 for quadratic LIV [166]. 
All these limits assumed s+ = 1. However, the Y?/dof was still very much greater 
than one when all the three datasets were combined. Various other variants were 
tried in this work (such as a constant lag in conjunction with intrinsic scatter, two 
different exponents in the intrinsic emission for GRB 190114C), but the reduced y? 
was still much greater than one for these other combinations also. 


16 Shantanu Desai 


Another search for LIV (for s+ = —1) using 8 Fermi-LAT GRBs with ener- 
gies above 100 MeV located at redshifts between 0.34 and 4.35 has also been re- 
cently carried out [167]. This work considered three distinct non-parametric mea- 
sures based on irregularity, skewness, and kurtosis to characterize the GRB emis- 
sion. A constant lag at the source frame was used for the intrinsic emission. This 
work obtained a 95% lower limit of Egg > 8.4 x 10!” GeV [167]. If one takes 
into account the intrinsic temporal spectral variations, the 95% lower limit becomes 
Egg > 2.4 x 10!" GeV [167]. 

A novel search for LIV using a forward modelling technique was carried out 
using the BATSE GRB catalog [168]. This work used the GRB catalog compiled 
in [169] after culling the GRBs with pseudo-redshifts greater than 4.5. They con- 
sidered two LIV theories. In the first one, the photon speed receives quadratic cor- 
rections due to the photon mass. In the second scenario, the photon speed depends 
linearly on the energy (similar to the n = 1 limit in Eq. 1.1). The expected time 
delays of photons were obtained using a novel forward modelling approach. This 
work obtained a limit on photon mass m, « 4.0 x 10754 eV and a quantum gravity 
length scale log < 5.3 x 1075 GeV-! at 95% c.l. (where log = Egg [168] These 
results are also agnostic to the astrophysical and observational uncertainties related 
to the spectral lag. These limits are also valid for both subluminal and superluminal 
propagation, where /og is negative. 

Most recently, a comprehensive search using 135 long GRBs from the 2020 
Fermi/GBM catalog was carried out [170]. Evidence for spectral lag transitions 
(from positive to negative) have been found for 32 of these GRBs in the redshift 
range 0.347-4.35 [170]. For each GRB, close to 20 spectral lags between 10 keV to 
10,000 keV (cf. Table 1 of [170]) were used to determine the transition. A broken 
power law was used to model the intrinsic spectral lag, which is given by: 


_ a - 1/u 
Aim =¢ (5#) des iO (528) | (1.11) 


where € denotes the normalization amplitude, o and o: represent the pre- and post- 
transition slopes before and after the transition energy Ep and u characterize the 
smoothness of the transition. Eq. 1.11 reduces to Eq. 1.10 when o = oo. Evidence 
for spectral lag transitions (from positive to negative) have been found using 32 
Fermi/GBM GRBs in the redshift range 0.347-4.35 [170]. The results from LIV 
search from each of the 32 GRBs for s+ = 1 can be found in Table 3 of [170]. 
The statistical significance of evidence for LIV based on the observed spectral lag 
transitions however has not been carried out for these GRBs. However, the authors 
have argued that some of the negative lags are astrophysical. For the linear LIV 
model, the lower limits for Egg range from (8.2 x 10!2-5.5 x 10? GeV [170]. For 
the quadratic LIV model, the corresponding limits for Egg range from (6.2 x 104 — 
1.7 x 107 GeV) [170]. Limits on the anisotropic LIV coefficients from SME were 
also set using this data [171]. 

Finally, observations of two potentially lensed GRBs (GRB 950830 and GRB 
200716C) have also been used to search for LIV [172]. Since lensing is achromatic, 
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the lensed images should have the same hardness and magnification. However, dif- 
ferent hardness values were observed in the energy bands, which could be induced 
due to LIV as argued in [172]. The LIV induced time delays due to lensing were 
then calculated and is energy dependent. For GRB 950830, the lens was estimated to 
have a redshift z ~ 1 and the source at redshift of about two. For GRB 200716C, the 
lens redshift was estimated to be z ~ 0.174 and the source at z ~ 0.348. The observed 
time delays were then used to constrain LIV assuming a point source lens mass. The 
best-fit values for a linear model of LIV and s+ = 1 is given by Egg > 3.2 x 10° 
GeV (GRB 950830) and Egg = 6.9 x 10° GeV (GRB 200716C) [172]. 


1.7.2 Searches using fixed source frame energies 


All the above studies have been done based on spectral lags in a fixed energy interval 
in the observer frame. However, if the data contains GRBs with different redshifts, 
the corresponding energy band in the source frame would be different. It has been 
shown that is a large scatter in the correlation between the observer-frame lags and 
the source-frame lags for the same GRB sample, implying that the observer-frame 
lags do not represent the rest-frame lag [173]. Therefore, the observer-frame lags 
would be strongly biased, since they sample different parts of the intrinsic light 
curves. Therefore, this could introduce systematics in the search for LIV [174, 20]. 
To ameliorate this, the search for LIV was done by calculating the spectral lags 
between two fixed energy intervals in the intrinsic or source frame [174]. Therefore, 


Eq. 1.5 and Eq. 1.9 get modified for a linear LIV and s+ = —1 as follows [174]: 
E! =F Zz 1 / 
eg mci f en (1.12) 
EgcHo(1+z) Jo A(z’) 


where Ej, and Ej are the rest frame energies in the upper and lower energy inter- 
vals, respectively. For this analysis, a dataset of 56 rest-frame spectral lags from 
the SWIFT satellite compiled in [175] was used [174]. This work also assumed 
a constant intrinsic lag (in the source frame). The spectral lags were again fit to a 
linear regression (similar to [2] and other aforementioned works) with the Lorentz 
violating term corresponding to the slope and the intrinsic delay corresponding to 
the intercept. The best-fit value for slope was consistent with 0 to within lo, im- 
plying there is no evidence for LIV. The 95% c.l. lower limit for Egg is given by 
Egg > (2 —2.2) x 10!^ GeV for s+ = —1 [174]. This was also the first work which 
divided the GRB sample into four groups according to redshift to look for a redshift- 
dependent trend in the LIV coefficients. No such variation was seen after splitting 
the data. An independent search for evidence for LIV using this dataset for s+ = 1 
and using the same intrinsic emission model as Eq. 1.10 was also carried out [166]. 
No evidence for either of the two LIV models was found using Bayesian model 
comparison tests [166]. 
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Most recently a search for LIV using a catalog of spectral lags from 44 short 
GRBs and 2 long GRBs from both SWIFT and FERMI-GBM was carried out [176]. 
The lags were constructed between the fixed source frame energy intervals of 15-70 
keV and 120-250 keV, and were obtained using a novel cross-correlation method, 
which utilizes the temporal information in the light curve and is agnostic to the 
details of the cross-correlation function [177, 178]. A constant source frame in- 
trinsic lag was used. Both AIC and BIC were used to assess the significance of a 
LIV signature. The intrinsic constant lag model as well as the LIV models provide 
equally good fits to the data. From this analysis, a lower limit on Egg in the range 
10!5 — 10!7 GeV was obtained at 95% c.l., depending on the choice of s [176]. 
A refined analysis has also been carried out using the energy-dependent intrinsic 
lag (Eq. 1.10), chronometers to characterize the expansion history, and also incor- 
porating the uncertainties in the energy intervals with these changes [179]. Unlike 
[176], no closed contours were obtained for Egg for both the linear and quadratic 
LIV hypotheses. The 95% c.l. lower limits obtained are Egg > 4 x 10!5 GeV and 
Egg > 6.8 x 10° GeV, respectively for s+ = 1 [179]. 

In addition to LIV with a pure sample of GRBs or a pure sample of AGN, a search 
for LIV using a combined sample of GRBs and AGNs has also been done using 
a combination of 24 GRBs and AGNs detected by Fermi-LAT, RHESSI, and TeV 
detectors (such as HESS, MAGIC and WHIPPLE) [180]. This data was used to con- 
strain 25 coefficients of the anisotropic birefringent Lorentz violating SME [180]. 


1.7.3 Searches based on TeV-PeV observations 


All the previous results for LIV were based on spectral lag measurements for GRBs 
in the keV-GeV energy regime (except for the MAGIC result which used TeV ob- 
servations of one GRB). Here, we discuss LIV searches based on one GRB detected 
at energies above 100 TeV, which has generated a lot of interest and excitement. 

In October 2022, GRB 2210094 was detected at energies at 18 TeV by the 
LHAASO observatory [181] and upto 250 TeV by the Carpet-2 detector [182]. 
This burst was also seen by Fermi-LAT in the 0.1-1.0 GeV range. This burst is 
also the brightest in terms of its peak flux and fluence [183]. The lower limit 
on the flux observed at LHAASO based on the detection of 1 photon is = > 
2.9 x 107? ph cm? s^! GeV~!. Extrapolating the Fermi-LAT spectrum to this en- 
ergy range gives a flux of 9.3 x 107!? ph cm~? s^! GeV~! . In a similar vein, the 
lower limit on the flux detected by Carpet-2 at 250 TeV is given by — 218x 
10716 ph cm? s^! GeV~!. The corresponding value extrapolated from Fermi-LAT 
is given by 6.7 x 107!4 ph cm? s^! GeV-!. Based on a comparison of the ob- 
served a interpolated spectrum at 18 TeV, a 95% c.l. upper limit on the opacity (7) , 
given by T < 17 was obtained [184]. This constraint is consistent with only a small 
fraction of models of extragalactic background light [185]. The corresponding limit 
at 250 TeV is given by T < 5.9, which is in severe tension with most models of 
extragalactic background light. This gives a upper limit on Egg, which is given by 
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Egg < 49Ep (for linear LIV) and Egg < 10-6E pl for quadratic LIV. Therefore, 
taken at face value, this also constitutes evidence for LIV. Other works which have 
analyzed this GRB have also come to similar conclusions [185, 186, 187, 188]. 
Again this result is in conflict with the lower limits obtained in [141]. However, this 
result also has a number of caveats associated with it as discussed in [184], from 
prosaic astrophysical assumptions [189] to other exotic Physics, which could ex- 
plain the disappearance of photons [190]. This result should soon be confirmed with 
more detections from LHAASO, HAWC, and the upcoming Cherenkov Telescope 
Array, which could detect VHE emission from GRBs upto 100s of TeV. 


1.8 Searches using polarization observations 


In addition to searches for LIV with photon counts, corresponding searches have 
also been carried out from polarization observations using GRBs as well as with 
AGNs. We now recap some of these results. 


1.8.1 Searches with GRBs 


The first test for LIV using polarization measurements was carried out for GRB 
0412194 located at a redshift of z = 0.02 [191]. In the presence of LIV, the pho- 
ton dispersion relation for frequency (@) and wavevector (k) gets modified to the 
following form: 

2E? 

Mpi ' 


o = K SE 


(1.13) 


where € is the LIV term and the rotation of direction of polarization after traversing 
a distance d is given by: 


kd 
AO(p) zs (1.14) 
r 
In terms of cosmological parameters this can be written as: 
H (+z) 
A0(p)= d ad 1.15 
(p) dag o A) (1.15) 


The two different values of o cause a rotation of the polarization during the prop- 
agation of linearly polarized photons. For this GRB a large amount of polarization 
was observed [192]. A search for a polarimetric angle shift between 200-250 keV 
and 250-325 keV, (which could be signature of LIV) was carried out. No such shift 
was seen. Based on this a constraint on the coupling & of the dimension five Lorentz 
violating interaction, which is given by é < 107^ [191]. 
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The subsequent search for LIV using polarized emission was done using GRB 
110721A, GRB 100826A and GRB 110301A using polarized emission from the 
IKAROS satellite [193]. The polarized emission was detected in two separate energy 
bands of 70-100 keV and 100-300 keV. A pseudo-redshift was estimated based on 
using the tightness of Epeak- Luminosity peak correlation. For constraints on LIV 
the 26 lower limit on the redshift was used. The most stringent limit was obtained 
for GRB 110721A with |E| < 2 x 107? — 8 x 10716, depending on the assumption 
used for the amount of intrinsic polarization [193]. 

Subsequently, a search for LIV using polarization measurements by ASTROSAT 
of 7 GRBs [194], as well as all of the five previously detected polarized GRBs was 
carried out [195] by applying a new technique proposed in [196]. For GRBs with 
unconfirmed redshifts, lower limits on each of the GRB redshifts were obtained 
by assuming the GRBs fall within 26 of the Amati relation [197, 198]. This work 
pointed out that all previous LIV-based limits based on polarized GRBs, which were 
obtained by assuming that the rotation angle of the polarization vector (A0), given 
by A0 < z/2 could be ameliorated by applying the technique proposed in [196]. 
Subsequently new limits were set using each of the 12 GRBs (cf. Table 1 in [195]) 
with the most stringent limit given by é < 5.2 x 10^" for GRB 061122 [195] 

In addition to polarization from the prompt emission phase of GRBs in the high 
energy range, searches for LIV from using polarization observations in the afterglow 
phase have also been carried out [199]. This work used optical afterglow observa- 
tions of GRB 020813 and GRB 021004 to look for birefringence. For GRB 020813, 
the best-fit value of € was € = 1.58 x 1077 and is consistent with a zero value to 
within 2.56 [199]. The corresponding value for GRB 021004 was € = ct at 
the 3o level [199]. 

More stringent limits on € can be obtained from upcoming missions such as 
POLAR, TSUBAME, COSI, and GRAPE which will carry out polarimetry of the 
prompt emission phase of the GRBs [200]. Other possible avenues for searching for 
LIV using astrophysical polarization measurements are reviews in [180]. 


1.8.2 Searches with AGNs 


A search for LIV-induced birefringence caused by rotation in the polarization vec- 
tor of a linearly polarized plane wave has also been done using optical spectropo- 
larimetry from five blazars (3C 66A, S5 0716+714, OJ 287, MK 421, and PKS 
2155-304) [201]. Polarization measurements for these aforementioned blazars were 
carried out using 37 groups of observations in five optical filters U BVRI. The 20 
limits on the five blazars are given by —8.91 x 1077 < é < 2.93 x 1075 [201]. 

A similar proof of principles test using observations of two AGNs (BL Lacartae 
and S5 B0716+714) was also carried out with the Array Photo polarimeter and 
constraints on Lorentz violation SME were set based on these observations [202]. 
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1.9 Conclusions 


In this manuscript, we have reviewed observational tests of LIV from neutrinos (us- 
ing terrestrial and astrophysical observations), gravitational waves, photon based 
astrophysical observations of GRBs, AGNs, and pulsars, and cosmological tests 
with galaxy clusters, BAO, Type la SNe and CMB. Most of the tests using GRBs 
and AGNs have been carried out using the observed spectral lags (in the keV-GeV 
energy range) by looking for an energy-dependent speed of light. Tests with po- 
larization measurements entail looking for a change in polarization angle due to 
vaccuum birefringence. Among these, some works have found evidence for such an 
energy dependence based on fitting the spectral lags, although not for the same value 
of s+ [64, 143, 145, 144, 148, 164, 160, 176]. Evidence for LIV based on TeV-PeV 
observations of GRB 221009A has also been claimed (although with caveats) [184]. 
However, the inferred values of Egg are in conflict some of the most stringent limits 
obtained [141, 142, 68] (See Table 3 of [21]). Furthermore, when one combines 
the data from disparate datasets, which individually point to LIV, one does not get 
a consistent picture or a good fit to the combined data [166]. Also, it is not easy to 
disentangle the LIV-induced lags from astrophysical contributions to the lag, which 
is still unknown. Many of the aforementioned works have used different models for 
the astrophysical lag. Given the diversity in the observed GRB light curves, one 
would need to incorporate these observational based priors, while modelling the in- 
trinsic spectral lags in order to make a robust case for evidence of LIV. The recent 
observation of GRB 221009A at energies above 100 TeV is exciting. More observa- 
tions of GRBs at such high energies should soon be forthcoming from high energy 
gamma ray telescopes at energies above > 100 GeV [203]. 

In the area of Cosmology the recent evidence for birefringence is tantalizing, and 
it would be interesting to see how it pans out with additional data from upcoming 
CMB missions [113]. 


Acknowledgements I would like to thank Cosimo Bambi and Alejandro Cardinas for invitation 
to write this chapter, and Shalini Ganguly, Haveesh Singirikonda and Rajdeep Agrawal for prior 
collaboration in this area. 


References 


[1] G. Amelino-Camelia, J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, S. Sarkar, 
Nature 395, 525 (1998). DOI 10.1038/26793 

[2] J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, A.S. Sakharov, E.K.G. Sark- 
isyan, Astroparticle Physics 25, 402 (2006). DOI 10.1016/j.astropartphys. 
2006.04.001 

[3] G. Amelino-Camelia, J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, S. Sarkar, 
Nature 393(6687), 763 (1998). DOI 10.1038/31647 


22 Shantanu Desai 


[4] A. Addazi, et al., Prog. Part. Nucl. Phys. 125, 103948 (2022). DOI 10.1016/ 
j-ppnp.2022.103948 
[5] J.P. Norris, R.J. Nemiroff, J.T. Bonnell, J.D. Scargle, C. Kouveliotou, W.S. 
Paciesas, C.A. Meegan, G.J. Fishman, Astrophys. J. 459, 393 (1996). DOI 
10.1086/176902 
[6] U. Jacob, T. Piran, JCAP 1, 031 (2008). DOI 10.1088/1475-7516/2008/01/ 
031 
[7] N. Aghanim, et al., Astron. Astrophys. 641, A6 (2020). DOI 10.1051/ 
0004-636 1/201833910 
[8] M. Biesiada, A. Piórkowska, Classical and Quantum Gravity 26(12), 125007 
(2009). DOI 10.1088/0264-9381/26/12/125007 
[9] G. Rosati, G. Amelino-Camelia, A. Marciano, M. Matassa, Phys. Rev. D. 
92(12), 124042 (2015). DOI 10.1103/PhysRevD.92.124042 
[10] V.A. Kostelecky, Phys. Rev. D. 69(10), 105009 (2004). DOI 10.1103/ 
PhysRevD.69.105009 
[11] V.A. Kostelecky, M. Mewes, Astrophys. J. Lett.689(1), L1 (2008). DOI 
10.1086/595815 
[12] S. Sarkar, Modern Physics Letters A 17(15-17), 1025 (2002). DOI 10.1142/ 
$0217732302007521 
[13] J. Bolmont, A. Jacholkowska, Advances in Space Research 47(2), 380 
(2011). DOI 10.1016/j.asr.2010.07.025 
[14] J. Ellis, N.E. Mavromatos, Astroparticle Physics 43, 50 (2013). DOI 10.1016/ 
j.astropartphys.2012.05.004 
[15] J. Bolmont, C. Perennes, J. Phys. Conf. Ser. 1586(1), 012033 (2020). DOI 
10.1088/1742-6596/1586/1/012033 
[16] D. Horns, A. Jacholkowska, Comptes Rendus Physique 17(6), 632 (2016). 
DOI 10.1016/j.crhy.2016.04.006 
[17] P. Petitjean, FY. Wang, X.F. Wu, J.J. Wei, Space Science Reviews202(1-4), 
195 (2016). DOT 10.1007/s11214-016-0235-6 
[18] E Stecker, Symmetry 9(10), 201 (2017). DOI 10.3390/sym9100201 
[19] H. Martinez-Huerta, R.G. Lang, V. de Souza, Symmetry 12(8), 1232 (2020). 
DOI 10.3390/sym1208 1232 
[20] J.J. Wei, X.F. Wu, Frontiers of Physics 16(4), 44300 (2021). DOI 10.1007/ 
s11467-021-1049-x 
[21] J.J. Wei, SE Wu, arXiv e-prints arXiv:2111.02029 (2021) 
[22] J.J. Wei, arXiv e-prints arXiv:2210.02114 (2022) 
[23] P. He, B.Q. Ma, Universe 8(6), 323 (2022). DOI 10.3390/universe8060323 
[24] C. Pérez de los Heros, T. Terzić, arXiv e-prints arXiv:2209.06531 (2022) 
[25] R.G. Lang, H. Martínez-Huerta, V. de Souza, Universe 8(8), 435 (2022). DOI 
10.3390/universe8080435 
[26] T. Jacobson, S. Liberati, D. Mattingly, Annals of Physics 321(1), 150 (2006). 
DOI 10.1016/j.aop.2005.06.004 
[27] J.D. Tasson, Reports on Progress in Physics 77(6), 062901 (2014). DOI 
10.1088/0034-4885/77/6/062901 
[28] J. Tasson, Symmetry 8(11), 111 (2016). DOI 10.3390/sym8110111 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 23 


[29] S. Liberati, Classical and Quantum Gravity 30(13), 133001 (2013). DOI 
10.1088/0264-938 1/30/13/133001 

[30] S. Liberati, L. Maccione, Annual Review of Nuclear and Particle Science 
59(1), 245 (2009). DOI 10.1146/annurev.nucl.010909.083640 

[31] D. Mattingly, Living Reviews in Relativity 8, 5 (2005). DOI 10.12942/ 
Irr-2005-5 

[32] G. Amelino-Camelia, Living Reviews in Relativity 16, 5 (2013). DOI 10. 
12942/1rr-2013-5 

[33] C.M. Will, in Einstein, 1905-2005: Poincar&eacute; Seminar 2005, vol. 47, 
ed. by T. Damour, O. Darrigol, B. Duplantier, V. Rivasseau (2006), p. 33. 
DOI 10.1007/3-7643-7436-5.2 

[34] B. Altschul, EPL (Europhysics Letters) 135(4), 41001 (2021). DOI 10.1209/ 
0295-5075/ac2ea2 

[35] P.D. Group, P. Zyla, R. Barnett, J. Beringer, O. Dahl, D. Dwyer, D. Groom, 
C.J. Lin, K. Lugovsky, E. Pianori, et al., Progress of Theoretical and Experi- 
mental Physics 2020(8), 083C01 (2020) 

[36] V.A. Kostelecky, J.D. Tasson, Physics Letters B 749, 551 (2015). DOI 10. 
1016/j.physletb.2015.08.060 

[37] S. Boran, S. Desai, E.O. Kahya, R.P. Woodard, Phys. Rev. D. 97(4), 041501 
(2018). DOI 10.1103/PhysRevD.97.041501 

[38] K. Bora, S. Desai, JCAP 2021(2), 012 (2021). DOI 10.1088/1475-7516/ 
2021/02/012 

[39] S. Bhagvati, S. Desai, Classical and Quantum Gravity 39(1), 017001 (2022). 
DOI 10.1088/1361-6382/ac3c8c 

[40] D.R. Lorimer, M. Kramer, Handbook of Pulsar Astronomy (2012) 

[41] T.T. Reddy Ch., S. Desai, New Astronomy 91, 101673 (2022). DOI 10.1016/ 
j.newast.2021.101673 

[42] R.D. Blandford, Philosophical Transactions of the Royal Society of London 
Series A 341(1660), 177 (1992). DOI 10.1098/rsta. 1992.0092 

[43] S. Desai, E.O. Kahya, Modern Physics Letters A 31(13), 1650083 (2016). 
DOI 10.1142/S02177323 16500838 

[44] P. Kaaret, Astron. & Astrophys. 345, L32 (1999) 

[45] M.L. Ahnen, et al., Astrophys. J. Suppl. 232(1), 9 (2017). DOI 10.3847/ 
1538-4365/aa8404 

[46] M. Amenomori, et al., Phys. Rev. Lett. 123(5), 051101 (2019). DOI 10.1103/ 
PhysRevLett.123.051101 

[47] A.U. Abeysekara, et al., Astrophys. J. 881, 134 (2019). DOI 10.3847/ 
1538-4357/ab2f7d 

[48] P. Satunin, European Physical Journal C 79(12), 1011 (2019). DOI 10.1140/ 
epjc/s10052-019-7520-y 

[49] A. Albert, et al., Phys. Rev. Lett. 124(13), 131101 (2020). DOI 10.1103/ 
PhysRevLett.124.131101 

[50] Z. Cao, et al., Science 373(6553), 425 (2021). DOI 10.1126/science.abg5137 

[51] C. Li, B.Q. Ma, Physics Letters B 829, 137034 (2022). DOI 10.1016/j. 
physletb.2022.137034 


24 Shantanu Desai 


[52] Z. Cao, et al., Nature 594(7861), 33 (2021). DOI 10.1038/ 
s41586-021-03498-z 

[53] Z. Cao, et aL, Phys. Rev. Lett. 128(5), 051102 (2022). DOI 10.1103/ 
PhysRevLett.128.051102 

[54] C. Li, B.Q. Ma, Phys. Rev. D. 104(6), 063012 (2021). DOI 10.1103/ 
PhysRevD.104.063012 

[55] EW. Stecker, S.L. Glashow, Astroparticle Physics 16(1), 97 (2001). DOI 
10.1016/S0927-6505(01)00137-2 

[56] T. Tanimori, K. Sakurazawa, S.A. Dazeley, P.G. Edwards, T. Hara, 
Y. Hayami, S. Kamei, T. Kifune, T. Konishi, Y. Matsubara, T. Matsuoka, 
Y. Mizumoto, A. Masaike, M. Mori, H. Muraishi, Y. Muraki, T. Naito, S. Oda, 
S. Ogio, T. Osaki, J.R. Patterson, M.D. Roberts, G.P. Rowell, A. Suzuki, 
R. Suzuki, T. Sako, T. Tamura, G.J. Thornton, R. Susukita, S. Yanagita, 
T. Yoshida, T. Yoshikoshi, Astrophys. J. Lett.492(1), L33 (1998). DOI 
10.1086/311077 

[57] J. Nishimura, M. Fujii, T. Taira, E. Aizu, H. Hiraiwa, T. Kobayashi, K. Niu, 
I. Ohta, R.L. Golden, T.A. Koss, Astrophys. J. 238, 394 (1980). DOI 10. 
1086/157997 

[58] C.M. Urry, P. Padovani, Pub. Astro. Soc. Pacific 107, 803 (1995). DOI 

10.1086/133630 

S.D. Biller, A.C. Breslin, J. Buckley, M. Catanese, M. Carson, D.A. Carter- 

Lewis, M.F. Cawley, D.J. Fegan, J.P. Finley, J.A. Gaidos, A.M. Hillas, 

F. Krennrich, R.C. Lamb, R. Lessard, C. Masterson, J.E. McEnery, B. McK- 

ernan, P. Moriarty, J. Quinn, H.J. Rose, F. Samuelson, G. Sembroski, P. Skel- 

ton, T.C. Weekes, Phys. Rev. Lett. 83(11), 2108 (1999). DOI 10.1103/ 

PhysRevLett.83.2108 

[60] J. Albert, et al., Phys. Lett. B 668, 253 (2008). DOI 10.1016/j.physletb.2008. 
08.053 

[61] F. Aharonian, et al., Phys. Rev. Lett. 101, 170402 (2008). DOI 10.1103/ 
PhysRevLett.101.170402 

[62] A. Abramowski, et al., Astropart. Phys. 34, 738 (2011). DOI 10.1016/j. 
astropartphys.2011.01.007 

[63] A. Abramowski, et al., Astrophys. J. 802(1), 65 (2015). DOI 10.1088/ 
0004-637X/802/1/65 

[64] L. Shao, Z. Xiao, B.Q. Ma, Astroparticle Physics 33(5-6), 312 (2010). DOI 
10.1016/j.astropartphys.2010.03.003 

[65] H. Abdalla, et al., Astrophys. J. 870(2), 93 (2019). DOI 10.3847/1538-4357/ 
aaflc4 

[66] H. Li, B.Q. Ma, Science Bulletin 65(4), 262 (2020). DOI 10.1016/j.scib. 
2019.11.024 

[67] J. Biteau, D.A. Williams, Astrophys. J. 812(1), 60 (2015). DOI 10.1088/ 
0004-637X/8 12/1/60 

[68] R.G. Lang, H. Martinez-Huerta, V. de Souza, Phys. Rev. D. 99(4), 043015 
(2019). DOI 10.1103/PhysRevD.99.043015 


[59 


— 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 25 


[69] J.S. Diaz, arXiv e-prints arXiv: 1406.6838 (2014). DOI 10.48550/arXiv. 1406. 
6838 

[70] J. Diaz, Symmetry 8(10), 105 (2016). DOI 10.3390/sym8 100105 

[71] Y. Fukuda, et al., Phys. Rev. Lett. 81, 1562 (1998). DOI 10.1103/ 
PhysRevLett.81.1562 

[72] Y. Ashie, et al., Phys. Rev. D 71, 112005 (2005). DOI 10.1103/PhysRevD. 
71.112005 

[73] S. Desai, et al., Phys. Rev. D 70, 083523 (2004). DOI 10.1103/PhysRevD. 
70.083523 

[74] S. Desai, et al., Astropart. Phys. 29, 42 (2008). DOI 10.1016/j.astropartphys. 
2007.11.005 

[75] S. Coleman, S.L. Glashow, Phys. Rev. D. 59(11), 116008 (1999). DOI 10. 
1103/PhysRevD.59.116008 

[76] G.L. Fogli, E. Lisi, A. Marrone, G. Scioscia, Phys. Rev. D. 60(5), 053006 
(1999). DOI 10.1103/PhysRevD.60.053006 

[77] K. Abe, et al., Phys. Rev. D 91(5), 052003 (2015). DOI 10.1103/PhysRevD. 
91.052003 

[78] D.F.G. Fiorillo, G. Mangano, S. Morisi, O. Pisanti, JCAP 04, 079 (2021). 
DOI 10.1088/1475-7516/2021/04/079 

[79] T. Adam, et al., JHEP 10, 093 (2012). DOI 10.1007/JHEP10(2012)093 

[80] M. Antonello, et al., Phys. Lett. B 713, 17 (2012). DOI 10.1016/j.physletb. 
2012.05.033 

[81] P. Adamson, et al., Phys. Rev. D 92(5), 052005 (2015). DOI 10.1103/ 
PhysRevD.92.052005 

[82] R. Ehrlich, Symmetry 14(6), 1198 (2022). DOI 10.3390/sym14061198 

[83] R.M. Bionta, et al., Phys. Rev. Lett. 58, 1494 (1987). DOI 10.1103/ 
PhysRevLett.58.1494 

[84] K. Hirata, et al., Phys. Rev. Lett. 58, 1490 (1987). DOI 10.1103/PhysRevLett. 
58.1490 

[85] M.J. Longo, Phys. Rev. D. 36(10), 3276 (1987). DOI 10.1103/PhysRevD.36. 
3276 

[86] R. Laha, Phys. Rev. D. 100(10), 103002 (2019). DOI 10.1103/PhysRevD. 
100.103002 

[87] S. Boran, S. Desai, E.O. Kahya, European Physical Journal C 79(3), 185 
(2019). DOI 10.1140/epjc/s10052-019-6695-6 

[88] J. Ellis, N.E. Mavromatos, A.S. Sakharov, E.K. Sarkisyan-Grinbaum, Physics 
Letters B 789, 352 (2019). DOI 10.1016/j.physletb.2018.11.062 

[89] K. Wang, S.Q. Xi, L. Shao, R.Y. Liu, Z. Li, Z.K. Zhang, Phys. Rev. D. 102(6), 
063027 (2020). DOI 10.1103/PhysRevD.102.063027 

[90] M.G. Aartsen, et al., Phys. Rev. Lett. 111, 021103 (2013). DOI 10.1103/ 
PhysRevLett.111.021103 

[91] E. Borriello, S. Chakraborty, A. Mirizzi, P.D. Serpico, Phys. Rev. D. 87(11), 
116009 (2013). DOI 10.1103/PhysRevD.87.116009 

[92] G. Amelino-Camelia, L. Barcaroli, G. D' Amico, N. Loret, G. Rosati, Physics 
Letters B 761, 318 (2016). DOI 10.1016/j.physletb.2016.07.075 


26 Shantanu Desai 


[93] G. Amelino-Camelia, G. D'Amico, G. Rosati, N. Loret, Nature Astronomy 
1, 0139 (2017). DOI 10.1038/s41550-017-0139 
[94] Y. Huang, H. Li, B.Q. Ma, Phys. Rev. D. 99(12), 123018 (2019). DOI 10. 
1103/PhysRevD.99.123018 
[95] Y. Huang, B.Q. Ma, Communications Physics 1(1), 62 (2018). DOI 10.1038/ 
s42005-018-0061-0 
[96] H. Zhang, L. Yang, Universe 8(5) 260 (2022). DOI 10.3390/ 
universe8050260 
[97] M.G. Aartsen, et al., Astrophys. J. 824(2), 115 (2016). DOI 10.3847/ 
0004-637X/824/2/115 
[98] M.G. Aartsen, et al., Astrophys. J. 843(2), 112 (2017). DOI 10.3847/ 
1538-4357/aa7569 
[99] B. Clark, JINST 16(10), C10007 (2021). DOI 10.1088/1748-0221/16/10/ 
C10007 
[100] M.G. Aartsen, et al., J. Phys. G 44(5), 054006 (2017). DOI 10.1088/ 
1361-6471/44/5/054006 
[101] L. Perivolaropoulos, F. Skara, New Astronomy Reviews 95, 101659 (2022). 
DOI 10.1016/j.newar.2022.101659 
[102] J.W. Moffat, International Journal of Modern Physics D 2(3), 351 (1993). 
DOI 10.1142/80218271893000246 
[103] J.D. Barrow, Phys. Rev. D. 59(4), 043515 (1999). DOI 10.1103/PhysRevD. 
59.043515 
[104] G. Rodrigues, C. Bengaly, JCAP 2022(7), 029 (2022). DOI 10.1088/ 
1475-7516/2022/07/029 
[105] LE.C.R. Mendonga, K. Bora, R.E.L. Holanda, S. Desai, S.H. Pereira, JCAP 
2021(11), 034 (2021). DOI 10.1088/1475-7516/2021/11/034 
[106] T. Liu, S. Cao, M. Biesiada, Y. Liu, Y. Lian, Y. Zhang, Mon. Not. R. Astron. 
Soc. 506(2), 2181 (2021). DOI 10.1093/mnras/stab1868 
[107] J. Magueijo, Astrophys. and Space Science 283(4), 493 (2003). DOI 10. 
1023/A:1022560802810 
[108] A. Lue, L. Wang, M. Kamionkowski, Phys. Rev. Lett. 83(8), 1506 (1999). 
DOI 10.1103/PhysRevLett.83.1506 
[109] E. Komatsu, et al., Astrophys. J. Suppl. 180, 330 (2009). DOI 10.1088/ 
0067-0049/180/2/330 
[110] T. Kahniashvili, R. Durrer, Y. Maravin, Phys. Rev. D. 78(12), 123009 (2008). 
DOI 10.1103/PhysRevD.78.123009 
[111] Y. Minami, E. Komatsu, Phys. Rev. Lett. 125(22), 221301 (2020). DOI 
10.1103/PhysRevLett.125.221301 
[112] P. Diego-Palazuelos, J.R. Eskilt, Y. Minami, M. Tristram, R.M. Sullivan, A.J. 
Banday, R.B. Barreiro, H.K. Eriksen, K.M. Górski, R. Keskitalo, E. Ko- 
matsu, E. Martínez-González, D. Scott, P. Vielva, I.K. Wehus, Phys. Rev. 
Lett. 128(9), 091302 (2022). DOI 10.1103/PhysRevLett.128.091302 
[113] E. Komatsu, Nature Reviews Physics 4(7), 452 (2022). DOI 10.1038/ 
s42254-022-00452-4 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 27 


[114] J.A. Zuntz, P.G. Ferreira, T.G. Zlosnik, Phys. Rev. Lett. 101(26), 261102 
(2008). DOI 10.1103/PhysRevLett.101.261102 

[115] L. Caloni, S. Giardiello, M. Lembo, M. Gerbino, G. Gubitosi, M. Lattanzi, 
L. Pagano, arXiv e-prints arXiv:2212.04867 (2022). DOI 10.48550/arXiv. 
2212.04867 

[116] B.P. Abbott, et al., Phys. Rev. Lett. 116(6), 061102 (2016). DOI 10.1103/ 
PhysRevLett.116.061102 

[117] D. Blas, M.M. Ivanov, I. Sawicki, S. Sibiryakov, Soviet Journal of Ex- 
perimental and Theoretical Physics Letters 103(10), 624 (2016). DOI 
10.1134/80021364016100040 

[118] B.P. Abbott, et al., Astrophys. J. Lett. 848(2), L13 (2017). DOI 10.3847/ 
2041-8213/aa920c 

[119] E.O. Kahya, S. Desai, Physics Letters B 756, 265 (2016). DOI 10.1016/. 
physletb.2016.03.033 

[120] P. Kumar, B. Zhang, Physics Reports 561, 1 (2015). DOI 10.1016/j.physrep. 
2014.09.008 

[121] A. Singh, S. Desai, JCAP 2022(2), 010 (2022). DOI 10.1088/1475-7516/ 
2022/02/010 

[122] C. Kouveliotou, C.A. Meegan, G.J. Fishman, N.P. Bhat, M.S. Briggs, T.M. 
Koshut, W.S. Paciesas, G.N. Pendleton, Astrophys. J. Lett.413, L101 (1993). 
DOI 10.1086/186969 

[123] S.E. Woosley, J.S. Bloom, Ann. Rev. Astron. Astrophys. 44, 507 (2006). DOI 
10.1146/annurev.astro.43.072103.150558 

[124] E. Nakar, Physics Reports 442, 166 (2007). DOI 10.1016/j.physrep.2007.02. 
005 

[125] S. Kulkarni, S. Desai, Astrophys. and Space Science 362, 70 (2017). DOI 
10.1007/s10509-017-3047-6 

[126] A. Bhave, S. Kulkarni, S. Desai, P.K. Srijith, Astrophys. and Space Science 
367(4), 39 (2022). DOI 10.1007/s10509-022-04068-z 

[127] B.E. Schaefer, Phys. Rev. Lett. 82(25), 4964 (1999). DOI 10.1103/ 
PhysRevLett.82.4964 

[128] J. Ellis, K. Farakos, N.E. Mavromatos, V.A. Mitsou, D.V. Nanopoulos, As- 
trophys. J. 535(1), 139 (2000). DOI 10.1086/308825 

[129] J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, A.S. Sakharov, Astron. & As- 
trophys. 402, 409 (2003). DOI 10.1051/0004-6361:20030263 

[130] S.E. Boggs, C.B. Wunderer, K. Hurley, W. Coburn, Astrophys. J. Lett.611(2), 
L77 (2004). DOI 10.1086/423933 

[131] J. Ellis, N.E. Mavromatos, D.V. Nanopoulos, A.S. Sakharov, E.K.G. Sark- 
isyan, arXiv e-prints arXiv:0712.2781 (2007). DOI 10.48550/arXiv.0712. 
2781 

[132] Y. Pan, Y. Gong, S. Cao, H. Gao, Z.H. Zhu, Astrophys. J. 808(1), 78 (2015). 
DOI 10.1088/0004-637X/808/1/78 

[133] X.B. Zou, H.K. Deng, Z.Y. Yin, H. Wei, Physics Letters B 776, 284 (2018). 
DOI 10.1016/j.physletb.2017.11.053 


28 Shantanu Desai 


[134] M. Rodriguez Martinez, T. Piran, Y. Oren, JCAP 2006(5), 017 (2006). DOI 
10.1088/1475-75 16/2006/05/017 

[135] R. Lamon, N. Produit, F. Steiner, General Relativity and Gravitation 40(8), 
1731 (2008). DOI 10.1007/s 107 14-007-0580-6 

[136] J. Bolmont, A. Jacholkowska, J.L. Atteia, F. Piron, G. Pizzichini, Astrophys. 
J. 676(1), 532 (2008). DOT 10.1086/527524 

[137] A.A. Abdo, et al., Science 323, 1688 (2009). DOI 10.1126/science.1169101 

[138] Z. Xiao, B.Q. Ma, Phys. Rev. D. 80(11), 116005 (2009). DOI 10.1103/ 
PhysRevD.80.116005 

[139] A.A. Abdo, M. Ackermann, M. Ajello, K. Asano, W.B. Atwood, M. Axels- 
son, L. Baldini, J. Ballet, G. Barbiellini, M.G. Baring, et al., Nature 462, 331 
(2009). DOI 10.1038/nature08574 

[140] R.J. Nemiroff, R. Connolly, J. Holmes, A.B. Kostinski, Phys. Rev. Lett. 
108(23), 231103 (2012). DOI 10.1103/PhysRevLett. 108.231103 

[141] V. Vasileiou, A. Jacholkowska, F. Piron, J. Bolmont, C. Couturier, J. Gra- 
not, F.W. Stecker, J. Cohen-Tanugi, F. Longo, Phys. Rev. D. 87(12), 122001 
(2013). DOI 10.1103/PhysRevD.87.122001 

[142] V. Vasileiou, J. Granot, T. Piran, G. Amelino-Camelia, Nature Physics 11, 
344 (2015). DOI 10.1038/nphys3270 

[143] S. Zhang, B.Q. Ma, Astroparticle Physics 61, 108 (2015). DOI 10.1016/. 
astropartphys.2014.04.008 

[144] H. Xu, B.Q. Ma, Astroparticle Physics 82, 72 (2016). DOI 10.1016/. 
astropartphys.2016.05.008 

[145] H. Xu, B.Q. Ma, Physics Letters B 760, 602 (2016). DOI 10.1016/j.physletb. 
2016.07.044 

[146] Y. Liu, B.Q. Ma, European Physical Journal C 78(10), 825 (2018). DOI 
10.1140/epjc/s10052-018-6294-y 

[147] M.G. Bernardini, G. Ghirlanda, S. Campana, P. D'Avanzo, J.L. Atteia, 
S. Covino, G. Ghisellini, A. Melandri, F. Piron, R. Salvaterra, G. Taglia- 
ferri, Astron. & Astrophys. 607, A121 (2017). DOI 10.1051/0004-6361/ 
201731332 

[148] J.J. Wei, B.B. Zhang, L. Shao, X.F. Wu, P. Mészáros, Astrophys. J. Lett.834, 
L13 (2017). DOI 10.3847/2041-8213/834/2/L13 

[149] L. Shao, B.B. Zhang, F.R. Wang, X.F. Wu, Y.H. Cheng, X. Zhang, B.Y. Yu, 
B.J. Xi, X. Wang, H.X. Feng, M. Zhang, D. Xu, Astrophys. J. 844(2), 126 
(2017). DOI 10.3847/1538-4357/aa7d01 

[150] C. Perennes, H. Sol, J. Bolmont, Astron. & Astrophys. 633, A143 (2020). 
DOI 10.105 1/0004-6361/201936430 

[151] S. Desai, EPL (Europhysics Letters) 115, 20006 (2016). DOI 10.1209/ 
0295-5075/115/20006 

[152] AR Liddle, Mon. Not. R. Astron. Soc. 351, L49 (2004). DOI 10.1111/j. 
1365-2966.2004.08033.x 

[153] A. Krishak, S. Desai, JCAP 2020(7), 006 (2020). DOI 10.1088/1475-7516/ 
2020/07/006 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 29 


[154] S. Ganguly, S. Desai, Astroparticle Physics 94, 17 (2017). DOI 10.1016/j. 
astropartphys.2017.07.003 

[155] G. Gunapati, A. Jain, P.K. Srijith, S. Desai, PASA 39, e001 (2022). DOI 
10.1017/pasa.2021.64 

[156] J.J. Wei, X.F. Wu, B.B. Zhang, L. Shao, P. Mészáros, V.A. Kostelecky, As- 
trophys. J. 842(2), 115 (2017). DOI 10.3847/1538-4357/aa7630 

[157] R.J. Lu, Y.F. Liang, D.B. Lin, J. Lü, X.G. Wang, H.J. Lü, H.B. Liu, E.W. 
Liang, B. Zhang, Astrophys. J. 865(2), 153 (2018). DOI 10.3847/1538-4357/ 
aadal6 

[158] S.S. Du, D.B. Lin, R.J. Lu, R.Q. Li, Y.Y. Gan, J. Ren, W. Xiang-Gao, EW 
Liang, Astrophys. J. 882(2), 115 (2019). DOI 10.3847/1538-4357/ab33fe 

[159] W.Q. Liang, R.J. Lu, C.F. Peng, W.H. Chen, Astrophys. J. 942(2), 67 (2023). 
DOI 10.3847/1538-4357/aca96d 

[160] Y. Pan, J. Qi, S. Cao, T. Liu, Y. Liu, S. Geng, Y. Lian, Z.H. Zhu, Astrophys. 
J. 890(2), 169 (2020). DOI 10.3847/1538-4357/ab6ef5 

[161] M. Seikel, C. Clarkson, M. Smith, JCAP 1206, 036 (2012). DOI 10.1088/ 
1475-75 16/2012/06/036 

[162] H. Singirikonda, S. Desai, European Physical Journal C 80(8), 694 (2020). 
DOI 10.1140/epjc/s10052-020-8289-8 

[163] R. Jimenez, A. Loeb, Astrophys. J. 573(1), 37 (2002). DOI 10.1086/340549 

[164] S.S. Du, L. Lan, J.J. Wei, Z.M. Zhou, H. Gao, L.Y. Jiang, B.B. Zhang, Z.K. 
Liu, X.F. Wu, E.W. Liang, Z.H. Zhu, Astrophys. J. 906(1), 8 (2021). DOI 
10.3847/1538-4357/abc624 

[165] V.A. Acciari, et al., Phys. Rev. Lett. 125(2), 021301 (2020). DOI 10.1103/ 
PhysRevLett.125.021301 

[166] R. Agrawal, H. Singirikonda, S. Desai, JCAP 2021(5), 029 (2021). DOI 
10.1088/1475-7516/2021/05/029 

[167] J. Ellis, R. Konoplich, N.E. Mavromatos, L. Nguyen, A.S. Sakharov, E.K. 
Sarkisyan-Grinbaum, Phys. Rev. D. 99(8), 083009 (2019). DOI 10.1103/ 
PhysRevD.99.083009 

[168] D.J. Bartlett, H. Desmond, P.G. Ferreira, J. Jasche, Phys. Rev. D. 104(10), 
103516 (2021). DOI 10.1103/PhysRevD.104.103516 

[169] H. Yu, S.Q. Xi, FY. Wang, Astrophys. J. 860(2), 173 (2018). DOI 10.3847/ 
1538-4357/aac2e3 

[170] Z.K. Liu, B.B. Zhang, Y.Z. Meng, Astrophys. J. 935(2), 79 (2022). DOI 
10.3847/1538-4357/ac81b9 

[171] J.N. Wei, Z.K. Liu, J.J. Wei, B.B. Zhang, X.F. Wu, Universe 8(10), 519 
(2022). DOI 10.3390/universe8100519 

[172] L. Lan, A. Piórkowska-Kurpas, X. Wen, M. Biesiada, K. Liao, H. Gao, Z. Li, 
Astrophys. J. 937(2), 62 (2022). DOI 10.3847/1538-4357/ac8df5 

[173] T.N. Ukwatta, K.S. Dhuga, M. Stamatikos, C.D. Dermer, T. Sakamoto, 
E. Sonbas, W.C. Parke, L.C. Maximon, J.T. Linnemann, P.N. Bhat, A. Eskan- 
darian, N. Gehrels, A.U. Abeysekara, K. Tollefson, J.P. Norris, Mon. Not. R. 
Astron. Soc. 419(1), 614 (2012). DOI 10.1111/j.1365-2966.2011.19723.x 


20 Shantanu Desai 


[174] J.J. Wei, SE Wu, Astrophys. J. 851(2), 127 (2017). DOI 10.3847/ 
1538-4357/aa9d8d 

[175] M.G. Bernardini, G. Ghirlanda, S. Campana, S. Covino, R. Salvaterra, J.L. 
Atteia, D. Burlon, G. Calderone, P. D’Avanzo, V. D’Elia, G. Ghisellini, 
V. Heussaff, D. Lazzati, A. Melandri, L. Nava, S.D. Vergani, G. Tagliaferri, 
Mon. Not. R. Astron. Soc. 446(2), 1129 (2015). DOI 10.1093/mnras/stu2153 

[176] S. Xiao, S.L. Xiong, Y. Wang, S.N. Zhang, H. Gao, Z. Zhang, C. Cai, Q.B. 
Yi, Y. Zhao, Y.L. Tuo, X.Q. Li, X.Y. Wen, Z.H. An, W.X. Peng, S.J. Zheng, 
F. Zhang, L.M. Song, T.P. Li, Astrophys. J. Lett.924(2), L29 (2022). DOI 
10.3847/2041-8213/ac478a 

[177] T.P. Li, J.L. Qu, H. Feng, L.M. Song, G.Q. Ding, L. Chen, ChJAA4, 583 
(2004). DOI 10.1088/1009-9271/4/6/583 

[178] S. Xiao, S.L. Xiong, S.N. Zhang, L.M. Song, EJ. Lu, Y. Huang, C. Cai, Q.B. 
Yi, X. Y. Song, W. Chen, M.Y. Ge, CZ. Liu, X.B. Li, C.K. Li, Y. Zhao, As- 
trophys. J. 920(1), 43 (2021). DOI 10.3847/1538-4357/ac1420 

[179] S. Desai, R. Agrawal, H. Singirikonda, European Physical Journal C 83(1), 
63 (2023). DOI 10.1140/epjc/s10052-023-11229-z 

[180] F. Kislat, H. Krawczynski, Phys. Rev. D. 92(4), 045016 (2015). DOI 10. 
1103/PhysRevD.92.045016 

[181] Y. Huang, S. Hu, S. Chen, M. Zha, C. Liu, Z. Yao, Z. Cao, T.L. Experiment, 
GRB Coordinates Network 32677, 1 (2022) 

[182] DD Dzhappuev, Y.Z. Afashokov, LM. Dzaparova, T.A. Dzhatdoev, 
E.A. Gorbacheva, I.S. Karpikov, M.M. Khadzhiev, N.F. Klimenko, A.U. 
Kudzhaev, A.N. Kurenya, A.S. Lidvansky, O.I. Mikhailova, V.B. Petkov, E.I. 
Podlesnyi, N.A. Pozdnukhov, V.S. Romanenko, G.I. Rubtsov, S.V. Troitsky, 
I.B. Unatlokov, I.A. Vaiman, A.F. Yanin, K.V. Zhuravleva, The Astronomer's 
Telegram 15669, 1 (2022) 

[183] E. Burns, et al., (2023) 

[184] J.D. Finke, S. Razzaque, Astrophys. J. Lett.942(1), L21 (2023). DOI 10. 
3847/2041-8213/acadel 

[185] A. Baktash, D. Horns, M. Meyer, arXiv e-prints arXiv:2210.07172 (2022). 
DOI 10.48550/arXiv.2210.07172 

[186] G. Zhang, B.Q. Ma, Chinese Physics Letters 40(1), 011401 (2023). DOI 
10.1088/0256-307X/40/1/011401 

[187] G. Galanti, M. Roncadelli, F. Tavecchio, arXiv e-prints arXiv:2210.05659 
(2022). DOI 10.48550/arXiv.2210.05659 

[188] H. Li, B.Q. Ma, arXiv e-prints arXiv:2210.06338 (2022). DOI 10.48550/ 
arXiv.2210.06338 

[189] S. Sahu, B. Medina-Carrillo, G. Sánchez-Colón, S. Rajpoot, Astrophys. J. 
Lett.942(2), L30 (2023). DOI 10.3847/2041-8213/acac2f 

[190] G. Galanti, M. Roncadelli, F. Tavecchio, arXiv e-prints arXiv:2211.06935 
(2022). DOI 10.48550/arXiv.2211.06935 

[191] P. Laurent, D. Gótz, P. Binétruy, S. Covino, A. Fernandez-Soto, Phys. Rev. 
D. 83(12), 121301 (2011). DOI 10.1103/PhysRevD.83.121301 


1 Astrophysical and Cosmological Searches for Lorentz Invariance Violation 31 


[192] S. McGlynn, D.J. Clark, A.J. Dean, L. Hanlon, S. McBreen, D.R. Willis, 
B. McBreen, A.J. Bird, S. Foley, Astron. & Astrophys. 466(3), 895 (2007). 
DOI 10.105 1/0004-6361:20066179 

[193] K. Toma, S. Mukohyama, D. Yonetoku, T. Murakami, S. Gunji, T. Mi- 
hara, Y. Morihara, T. Sakashita, T. Takahashi, Y. Wakashima, H. Yonemochi, 
N. Toukairin, Phys. Rev. Lett. 109(24), 241104 (2012). DOI 10.1103/ 
PhysRevLett.109.241104 

[194] T. Chattopadhyay, S.V. Vadawale, E. Aarthy, N.P.S. Mithun, V. Chand, 
A. Ratheesh, R. Basak, A.R. Rao, V. Bhalerao, S. Mate, B. Arvind, 
V. Sharma, D. Bhattacharya, Astrophys. J. 884(2), 123 (2019). DOI 
10.3847/1538-4357/ab40b7 

[195] J.J. Wei, Mon. Not. R. Astron. Soc. 485(2), 2401 (2019). DOI 10.1093/ 
mnras/stz594 

[196] H.N. Lin, X. Li, Z. Chang, Mon. Not. R. Astron. Soc. 463(1), 375 (2016). 
DOI 10.1093/mnras/stw2007 

[197] L. Amati, F. Frontera, M. Tavani, J.J.M. in't Zand, A. Antonelli, E. Costa, 
M. Feroci, C. Guidorzi, J. Heise, N. Masetti, E. Montanari, L. Nicastro, 
E. Palazzi, E. Pian, L. Piro, P. Soffitta, Astron. & Astrophys. 390, 81 (2002). 
DOI 10.1051/0004-6361:20020722 

[198] G. Govindaraj, S. Desai, JCAP 2022(10), 069 (2022). DOI 10.1088/ 
1475-7516/2022/10/069 

[199] J.J. Wei, X.F. Wu, European Physical Journal Plus 135(6), 527 (2020). DOI 
10.1140/epjp/s13360-020-00554-x 

[200] M.L. McConnell, New Astronomy Reviews 76, 1 (2017). DOI 10.1016/. 
newar.2016.11.001 

[201] Q.Q. Zhou, S.X. Yi, J.J. Wei, SE Wu, Galaxies 9(2), 44 (2021). DOI 10. 
3390/galaxies9020044 

[202] A.S. Friedman, D. Leon, K.D. Crowley, D. Johnson, G. Teply, D. Tytler, 
B.G. Keating, G.M. Cole, Phys. Rev. D. 99(3), 035045 (2019). DOI 10.1103/ 
PhysRevD.99.035045 

[203] B. Zhang, Nature 575(7783), 448 (2019). DOI 10.1038/d41586-019-03503-6 


